Quantum Brownian Motion in a Bath of Parametric 
Oscillators: A model for system-field interactions 



B. L. Hu* 

Department of Physics, University of Maryland, College Park, MD 20742, USA 

A. Mataczt 

Department of Physics, University of Adelaide, 5005, Australia 



o ■ 

op . 3 December 1993 

Q: 

m ■ (umdpp 93-210) 

^ ! 

Abstract 

J> ' The quantum Brownian motion paradigm provides a unified framework where one 

■ can see the interconnection of some basic quantum statistical processes like decoher- 
Q ! ence, dissipation, particle creation, noise and fluctuation. The present paper continues 
CN ' the investigation into these issues begun in two earlier papers by Hu, Paz and Zhang 

■ on the quantum Brownian motion in a general environment via the influence func- 
, tional formalism. Here, the Brownian particle is coupled linearly to a bath of the 
' most general time dependent quadratic oscillators. This bath of parametric oscillators 

mimics a scalar field, while the motion of the Brownian particle modeled by a single 
5—1 . oscillator could be used to depict the behavior of a particle detector, a quantum field 

mode or the scale factor of the universe. An important result of this paper is the 
^ ^ ' derivation of the influence functional encompassing the noise and dissipation kernels in 

^ . terms of the Bogolubov coefficients, thus setting the stage for the influence functional 

^ I formalism treatment of problems in quantum field theory in curved spacetime. This 

- - ' method enables one to trace the source of statistical processes like decoherence and dis- 

sipation to vacuum fiuctuations and particle creation, and in turn impart a statistical 
mechanical interpretation of quantum field processes. With this result we discuss the 
statistical mechanical origin of quantum noise and thermal radiance from black holes 
and from uniformly-accelerated observers in Minkowski space as well as from the de 
Sitter universe discovered by Hawking, Unruh and Gibbons-Hawking. We also derive 
the exact evolution operator and master equation for the reduced density matrix of 
the system interacting with a parametric oscillator bath in an initial squeezed thermal 
state. These results are useful for decoherence and backreaction studies for systems 
and processes of interest in semiclassical cosmology and gravity. Our model and results 
are also expected to be useful for related problems in quantum optics. 

*Email: hu@umdhep.umd.edu 
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1 Introduction 



In two earlier papers, called Paper I, II henceforth Paz, Zhang and one of us began 

a systematic study of the celebrated problem of quantum Brownian motion (QBM) in a 
general environment using the Feynman- Vernon influence functional (IF) formalism p, ^, 
The special features associated with a nonohmic bath, or ohmic bath at low temperatures 
are the appearance of colored noise and nonlocal dissipation. The motivation for this study 
was amply explained there. What prompted them to this undertaking was the belief that 
a correct and deepened understanding of many interesting quantum statistical processes 
in the early universe and black holes requires an extension of the existing framework 
of quantum field theory in curved spacetime to statistical and stochastic fields in the 
framework of quantum open systems represented by the QBM 0. This viewpoint and 
methodology have indeed been applied to the analysis of some basic issues in quantum 
cosmology |10|, |lT], [l^, |13|, |T5[ , effective field theory [|16], , and the foundation of quantum 
mechanics, such as the uncertainty principle [|l^, |l^] and, most significantly, decoherence 
n|, ^ in the quantum to classical transition problem. (See the recent reviews 



of m 



and references therein and in Papers I, II for the standard literature on this 
topic). QBM is one of the two major paradigms of non-equilibrium statistical mechanics (the 
other being Boltzmann's kinetic theory) which is also amenable to detailed analysis. The 
study of many problems mentioned above which have nonlinear and nonlocal characteristics 
typical of quantum processes in gravitation and cosmology necessitates a closer scrutiny of 
this model beyond the ordinary limited conditions. 

As stated in the Introductions of Papers I and II, in order to make it useful for addressing 
issues in semiclassical gravity and quantum cosmology, a theory of quantum open systems 
has to be developed for quantum fields in curved spacetime. Noticeable effort has been put 
into this direction. Hu, Paz and Zhang constructed a stochastic field theory based on 
the QBM model and described how thermal field theory can be obtained as the equilibrium 
limit. As a tool for the study of the quantum origin of noise, fluctuations and structure 
formation in cosmology, they have extended the result of Paper II to quantum fields in 
Minkowski, Robertson- Walker and de Sitter spacetime. The nature and origin of quantum 
noise from particle-field interaction were discussed in where a statistical field-theoretical 
derivation of thermal radiance in the Hawking pi], B2| and Unruh effects 



were given. 

For semiclassical gravity Kuo and Ford have studied the fluctuations of quantum fields 
on the Einstein equations. Calzetta and Hu and the present authors ||3^ have analyzed 
the nature of noise, fluctuations, particle creation and backreaction for quantum fields in 
cosmological spacetimes and proposed an Einstein- Langevin equation as the centerpiece of a 
generalized theory of semiclassical gravity. For quantum cosmology, Sinha and Hu had 
used the coarse-grained 



J8| Schwinger-Keldysh effective action |^ to analyze the validity of 
the minisuperspace approximation in quantum cosmology. Paz and Sinha had used the 
influence functional method to discuss the transition from quantum to semiclassical gravity, 
and Calzetta and Hu |^ have studied dissipation problem in quantum cosmology. However, 
except for the few cases mentioned above, none of these earlier work made use of the master 
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or Langevin equation approach characteristic of the QBM study, which is necessary to probe 
fluctuation [05|, I 



instabihty and phase transition jT^] aspects of quantum 



into the noise, 
fields and spacetime. 

The present paper is an intermediate step in that direction. It is a generalization of 
Papers I and II in that the oscillators which make up the system and bath are now the 
most general time- dependent quadratic oscillators. This bath of parametric oscillators (as 
the number of modes goes to infinity) is identical to a scalar field, while the motion of 
the Brownian particle modeled by a single oscillator could be used to depict the behavior 
of a particle detector (with zero spring constant, as in e.g., |^), the scale factor of the 
universe, (with a negative kinetic energy term, as is seen in Eq.(2.2) of or the ho- 
mogeneous or inhomogeneous (density fluctuation) modes of the infiaton field in an early 

m 



universe 



12|, 0, ^ ^ Indeed the results obtained here can be taken over 
directly for the description of scalar fields in cosmological spacetimes, as our examples will 
demonstrate. Parametric amplification of the bath oscillator quanta gives rise to particle 
creation, as was pointed out by Parker and Zel'dovich EHI, which can be depicted by the Bo- 



golubov transformation between the creation and annihilation operators of the Fock spaces 
defined at different times. The averaged effect of the bath on the system is described by the 
influence functional, which, in the statistical field-theory context measures the backreaction 
of quantum processes associated with the field like particle creation on the dynamics of the 
background spacetime [^, There are two components in the influence functional, a noise 
kernel and a dissipation kernel. The noise kernel governs the decoherence process and also 
limits the degree of attainment of classicality [^]. It also depicts the effect of fluctuations 
(in particle number) BH]. The dissipation kernel which appears in the effective equation 



of motion depicts the effect of particle creation on the dynamics of the system. The QBM 
paradigm thus provides a unified framework where one can see the interconnection of the 
basic quantum statistical processes like decoherence, dissipation, particle creation, noise and 
fluctuation. The neccesity of analyzing these processes on the same footing was emphasized 
earlier in |TT[] . 

An important result of this paper is the derivation of the influence functional and thus the 
noise and dissipation kernels in terms of the Bogolubov coefficients. This enables one to trace 
the source of statistical processes like decoherence and dissipation to vacuum fluctuations 
and particle creation, and in turn impart a statistical mechanical interpretation of quantum 
field processes. With this we discuss the well-known results by Unruh [p3|. Hawking pl| 



and Gibbons- Hawking on thermal radiance from uniformly-accelerated observers [^, 
black holes and for comoving observers in de Sitter spacetime. From the explicit form of 
the noise and dissipation kernels we derived, one can see clearly the interplay of different 
factors which contribute to making the spectrum of particle creation in these situations 
thermal, and, more interestingly, what makes them nonthermal, as in the more general 
non-equilibrium situations. This is where the capability of the statistical field-theoretical 
interpretation supercedes the geometric interpretation (in terms of event horizons). We will 
discuss the implications of this point later. 

Although we have used examples from quantum and semiclassical cosmology to illustrate 
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the physical relevance of the QBM model with parametric bath, the range of applicability of 
this model goes beyond. An important area where parametric amplification plays a central 
role is in quantum optics. Here the properties of baths prepared in squeezed initial states 
(rigged reservoirs) are of interest ||5^, Squeezed baths are capable of processing optical 
signals (attenuation or amplification) while retaining their quantum features. It has also been 
shown that an appropriately squeezed bath is capable of greatly increasing the decoherence 
timescale [Q. The description of these processes is based on the quantum optical master 



equation generalised to include squeezing in the initial state. It is an approximate equation 
derived under the rotating wave, Born and Markov approximations. Since our formalism 
is exact it is capable of a more accurate description of non-equilibrium quantum statistical 
processes in quantum optics. It also allows for the squeezing to be generated dynamically 
rather than imposed as an initial condition. 

The effect of the bath on the system is studied here, as in the previous two QBM papers, 
by means of the influence functional formalism. We will derive exactly the evolution operator 
for the reduced density matrix, the influence functional, and the master equation for a time- 
dependent system and bath, using a slightly different method and language from Paper I. We 



adopt the language of squeeze and rotation operators |55, 56, 57 for describing the evolution 



of the system. In Sec. 2 we define the model and mention its relevance in problems in 
quantum optics, quantum and semiclassical cosmology, and quantum field theory in curved 
spacetimes. We then derive an analytic expression for the influence functional of a system 
linearly coupled to a bath of parametric oscillators in terms of the Bogolubov coefficients. 
In Sec. 3 we derive the exact evolution operator for the reduced density matrix and adopt 
the simpler method introduced by Paz and used in for the derivation of the master 



equation. We consider the general case when the bath is initially in a squeezed thermal state, 
which includes the common cases of a thermal state and a squeezed vacuum. We indicate 
how it is different from the model with a bath of time independent oscillators. The diffusion 
coefficients of this equation can be analyzed for decoherence studies, as is done in Papers 
I, Refs. ||2^ and ^0[. The relation of decoherence and particle creation was also discussed 
in the field theory context by Calzetta and Mazzitelli |Q and in the quantum cosmology 
context by Paz and Sinha |T^. Here we aim not at the decoherence or the dissipation 
processes, but focus on the definition and nature of noise associated with quantum fields and 
use them to depict some well-known processes such as the Hawking effect in gravitation and 
cosmology. 

In Sec. 4 we give a few simple examples of a system interacting with a bath of parametric 
oscillator, first treating the case with constant frequency, but with an initial squeezed thermal 
state and then the case of inverted oscillators which can be used to model amplifiers in 
quantum optics and electronics |5^. In Sec. 5, we derive the noise kernels for four cases: the 
accelerated observer, a two-dimensional black hole and a massless, conformally and minimally 
coupled scalar field in the de Sitter universe. In the de Sitter universe case the parametric 
oscillator bath can serve as a relatively simple model of the environment for homogenous 
and inhomogenous (density fluctuation) modes of the inflaton field in the early universe. We 
show the factors entering into the determination of the spectrum, and indicate how one can 
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understand the Hawking and Unruh effects in a purely statistical-mechanical sense without 
recourse to geometric notions (like the event horizon). We will discuss the ffuctuation- 
dissipation relation approach [^1], ^ 33] to understanding backreaction in semiclassical 
gravity in later work [llf 



I 36 



In Sec. 6 we summarize our results and suggest further 
problems in cosmology and gravitation where our results can be usefully applied. The details 
of derivation in Sec. 2 are recorded in the Appendices. 



2 Influence Functional 

Our system, the Brownian particle, is modeled by a parametric oscillator with mass M{s), 
cross term B{s) and natural (bare) frequency Q{s). The environment (bath) is also modeled 
by a set of parametric oscillators with mass m„(s), cross term bn{s) and natural frequency 
u!n{s)- The system is coupled to the bath through an arbitrary function F{x) of the system 
variable and linear in the bath variables g„ with coupling strength c„(s) in each oscillator. 
The action of the combined system + environment is 



S[x, q] = S'[x] + Ssiq] + Sint[x, q] 

r ri 

= J ds -M{s)(^x'^ + B{s)xx-n\s] 



x^ 



(2.1) 



+ Y.{l^n{s){qn + bn{s)qnqn " ^^(s)?^)} + J] (-c„(s)F(x)g„) 



where x and g„ are the coordinates of the particle and the oscillators. The bare frequency 
fl is different from the physical frequency flp due to its interaction with the bath, which 
depends on the cutoff frequency. We will discuss this point in more detail in Sec. 4.1. For 
problems discussed here we are interested in how the environment affects the system in some 
averaged way. The quantity containing this information is the reduced density matrix of the 
system obtained from the full density operator of the system + environment by tracing out 
the environmental degrees of freedom. The evolution operator is responsible for the time 
evolution of the reduced density matrix. The equation of motion governing this reduced 
density matrix is the master equation. Our central task is to derive the evolution operator 
and the master equation for the Brownian particle in a general environment. 

We will briefly review here the Feynman- Vernon influence functional method for deriving 
the evolution operator. Readers who are familiar with it can skip this subsection. The 
method provides an easy way to obtain a functional representation for the evolution operator 
for the reduced density matrix jTr- Let us start first with the evolution operator for the full 
density matrix J' defined by 

p{t)=J{t,U)p{ti). (2.2) 
As the full density matrix p evolves unitarily under the action of (2.1), the evolution 
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operator J has a simple path integral representation. In the position basis, the matrix 
elements of the evolution operator are given by 



J{x, q, a; , q , t | Xi, q^, x^, q^, ti) = /C(x, q, t \ Xi, q^, ti)lC{x , q , t | x^, q^, U) 

X q . x' q! 

= J Dx J Dq exp^S[x,q] J Dx' J Dq' exp -^^[x', q'] 



(2.3) 



where the operator /C is the evolution operator for the wave functions. In the second equation, 
the path integrals are over all histories compatible with the boundary conditions. We have 
used q to represent the full set of environmental coordinates g„ and the subscript i to denote 
the initial variables. 

The reduced density matrix is defined as 



-oo +00 



Prix,x') = J dq y dg'p(a;,q;x',q')5(q-q') (2.4) 



-00 —00 



and is propagated in time by the evolution operator J'r 

+00 +00 

P'P^Xt^ X , i) — / dx ^ I dx ».^^j'(*^, X , t I Xj^t^ '^i' ^2) Pv^'^i'} jj'i )* ^^*^) 



-00 —00 



By using the functional representation of the full density matrix evolution operator given 
in (2.3), we can also represent JT^ in path integral form. In general, the expression is very 
complicated since the evolution operator Jr depends on the initial state. If we assume that 
at a given time t = ti the system and the environment are uncorrelated 

p{t = ti) = p,{ti) X pb{ti), (2.6) 

then the evolution operator for the reduced density matrix does not depend on the initial 
state of the system and can be written as 

J'r{xf,x'j,t I Xi,x'^,ti) = J Dx J Dx' exp^|S'[a;] — S'[x']| J-'[x,x'] 

x'- 

(2.7) 

= J Dx J Dx' exp —A[x, x'] 



X. 



where the subscript / denotes final variables, and A[x, x'] is the effective action for the open 
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quantum system. The factor jF[x,x'], called the 'influence functional', is defined as 



hco +00 +00 1/ q/ 



T[x,x'] = J dqf J dm J dc^. jDqJ Dc( 



-00 —00 —00 qi q'. 

i 

n 



X exp Tf'S'blq] + Sint[x, q] - Sb[q] - Sint[x' , q']}pfe(qi, q-, U 



{21 



= exp j^5A[x, x'] 

where is the influence action. Thus = S[x] — S[x'] + 

It is not difficult to show that (2.8) has the representation independent form 

T[x,x'] = Tr[U{t,U)pt{U)U'\t,U)] (2.9) 

where f/(t) and U'(t) are the quantum propagators for the actions 

SE[(i\ + Sint[x{s), q] and S E[(i\ + Sint[x' (s) , q] and x{s) and x'{s) are treated as time dependent 
classical forcing terms. We have found this form to be much more convenient for deriving 
the influence functional. 

It is obvious from its definition that if the interaction term is zero, the influence functional 
is equal to unity and the influence action is zero. In general, the influence functional is a 
highly non-local object. Not only does it depend on the time history, but -and this is 
the more important property- it also irreducibly mixes the two sets of histories in the path 
integral of (2.7). Note that the histories x and x' could be interpreted as moving forward and 
backward in time respectively. Viewed in this way, one can see the similarity of the influence 
functional and the generating functional in the closed-time-path (CTP or Schwinger- 



Keldysh) integral formalism The Feynman rules derived in the CTP method are very 
useful for computing the IF. We shall treat the field theoretic problems in later papers. 

In those cases where the initial decoupling condition (2.6) is satisfied, the influence func- 
tional depends only on the initial state of the environment. The influence functional method 
can be extended to more general conditions, such as thermal equilibrium between the system 



and the environment [^, or correlated initial states p|, 



We now proceed to derive the influence functional for the model (2.1). From its definition 
it is clear that the influence functional is independent on the choice of system but only on the 
coupling of the system to the environment. Since our method is quite general we have been 
able to include, in Appendix A, the influence functional for the most general coupling linear 
in the bath variable. However in the body of the paper we only consider the position-position 
coupling in (2.1). For the case of a squeezed thermal initial state (to be defined later) we 
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find that for tlie model (2.1) tlie influence functional has the form 



x] = exp 



ds J ds'[F{x{s)) - F{x'{s))\fi{s, s') [F{x{s')) + F{x'{s')) 

t s 

>(x(s)) - F{x'{s))] u{s, s') \f{x{s')) - F{x'{s')) 



(2.10) 



— — I ds I ds' 



ti 



The functions /i(s, s') and i/(s, s') contain the effects of the environment on the system. 
They are known respectively as the dissipation and noise kernels. The reason for these 
names becomes clear in the semi-classical regime of the open system generated by (2.10). 

To find the appropriate semiclassical limit of this open quantum system we must find an 
action which generates the same influence functional as (2.10). Consider the action 



S[a{s)] = J' ds(L{x, X, s) + F{x)C{s 



(2.11) 



where ^(s) is a gaussian stochastic force with a non-zero mean. This system generates the 
influence functional . ^ 

J^[S, A] = ^exp ^ ' ^{s)A{s)ds ^ (2.12) 

where S and A are given by 

1 



S(.) = -{Fix{s)) + F{x\s))), A(s) = F{x{s)) - F{x'{s)) 



(2.13) 



and the average is understood as a functional integral over ^(s) multiplied by a normalised 
gaussian probability density functional V[^{s)] S(s)]. The probablility density functional is a 
functional of E(s) if we allow the statistical properties of ^ to depend on the system history. 
The avergaing can be performed to give ||3^ 



. t t s 

T\x,x'\ = expj^ J dsA{s){^{s)) - ^ J ds J ds'A{s)A{s')a 



2{S,S 



(2.14) 



where 6*2(5, s') is the second cumulant of the force ^. The equation of motion generated by 
the action (2.11) is 



dL d dL ^ dFjx) 
dx dt dx dx 



dx 



(2.15) 



where ^(t) is a zero mean gaussian stochastic force with {^(t)^(t')) = 6*2(5,5'). Now by 
comparing (2.14) and (2.10) we see that 



(^(s)) = -2 / ds'^{s,s')T,{s'), C2{s,s') = hu{s,s'). 



(2.16) 
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Therefore the semiclassical equation for the system described by the influence functional 
(2.10) is 

where {i{t)^{t')) = hu^t, t'). Under special circumstances fi tends to the derivative of a delta 
function which generates local dissipation. More generally we see that in the semiclassical 
limit IX generates non-local dissipation while hv is the correlator of a zero mean gaussian 
stochastic force. 

We find that the dissipation and noise kernels take the form 



1 r 

li{s, s') = - dujl{uj, s, s') {a^{s) + (3Us)VMs') + PM')} 

2 Jo I 



(2.18) 



1 r°° 

u(s, s') = - / dujl(uj, s, s') coth 

2 Jo 



' huj{ti) 
2kBT 



X 



(2.19) 



cosh2r(a;){a,(s) + (3^{s)}* {a^{s') + (3^{s')} 

+ cosh2r{u;){a^{s)+(3^{s)}{a^{s') + p^{s')y 

- sinh2r(a;)e-^^'^(-){a,(s) + /3^(s)}*{a^(s') + P^{s')} 

- sinh2r(a;)e2^'^(-H«c.(5) + PM)}Ms') + PM')} 

The quantities in these kernels describe three different properties of the environment. 

A) The a and (3, known as the Bogolubov coefficients, are complex numbers that contain 
all the information about the quantum dynamics of the bath parametric oscillators. They 
are derived from two coupled first order equations 



where the time dependent coefficients are given by 

^ 1 / mn{t)ul(t) ^ mn(t)hl(t) _ 

2 \ Kn 4:Kn m„(t) 



(2.20) 



(2.21) 



These equations are a by product of finding the quantum propagator for a parametric os- 
cillator which is done in appendix B. We will usually choose k„ (defined by (A. 6)) so that 
f{ti) = 0. Thus if 6„ = we will usually have Hn = ^n{ti)uJn{ti) ■ Eq's (2.21) must satisfy 
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the initial conditions a{ti) = l,(3{ti) = 0. Note that the mode label u in the kernels is 
equivalent to n in the continuous limit. 

If we assume 6 = and m = 1 we can show using (2.20) that 



Xn + u;l{t)Xn = 

where Xn{t) = a„(t) + Pn{t). The solution of (2.22) must satisfy X„(tj 
the noise and dissipation kernels become 



(2.22) 
1. In this case 



iy[s, s 



1 

2 Jo 



/^(■s, ■s') = o / dujl{uj, s, s') 

I JO 



x:is)x4s') - x^{s)x:{s') 



(2.23) 



dujl{uj, s, s') coth 



2kBT 



cosh 2r(uj] 



x:is)xM') + xus)x:{s' 



- sinh2r(a;) [e-^'^'-^^ X:{s)X:{s') + e^'^'^'^^X^{s)X^{s') 
Note that we can always write 

Xr,{t)=Cn{t)-tUJ„iti)Sn{t) 



(2.24) 



(2.25) 



where C„ and 5"^ are subject to the boundary conditions Cniti) = S'„(tj) = 1 and S'„(tj) = 
Cn{ti) = 0. If the kernels are written in this notation we can show that for a thermal ini- 



tial state (2.17) reduces to the classical Langevin equation in the high temperature limit ^9 
B) The spectral density, I{u, s, s') defined formally by 



I{U, S, S') =^6{lJ - Un) 



Cn[S)Cn[S 



(2.26) 



is obtained in the continuum limit. It contains information about the environmental mode 
density and coupling strength as a function of frequency. Different environments are clas- 
sified according to the functional form of the spectral density I{uj). On physical grounds, 
one expects the spectral density to go to zero for very high frequencies. Let us introduce a 
certain cutoff frequency A (a property of the environment) such that I{uj) — for c<j > A. 
The environment is classified as ohmic |Q, |^ if in the physical range of frequencies {u < A) 
the spectral density is such that I{uj) ~ u, as supra-ohmic if I{uj) ~ tt;"',n > 1 or as 
sub-ohmic if n < 1. The most studied ohmic case corresponds to an environment which in- 
duces a dissipative force linear in the velocity of the system. We will show this in section 4.1. 



C) The initial state of the bath is a squeezed thermal state. It has the form 



(2.27) 
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where pth is a thermal density matrix of temperature T defined by (A. 18) and S{r,(f)) is a 
squeeze operator defined by (B.12). Since a squeezed thermal state is still gaussian it is a 
tractable generalisation of the usual thermal initial state that is of interest in quantum optics 
5^ . For the case of zero temperature we have a squeezed vacuum initial state. 

Physically the term squeezing arises because the phase space noise distribution of a 
squeezed vacuum is an eUipse squeezed from a circle (coherent state) by r and rotated by 
angle (j) with respect to the x and p axes. Thus, for the squeezed vacuum |^| 



(g^) = — [cosh 2r — sinh 2r cos 2(f)] 

'^^ (2.28) 
(p^) = — [cosh 2r + sinh 2r cos 20] . 

Note that the dissipation kernel is independent of the bath initial state. 

For the case of no initial squeezing we see that the noise and dissipation kernels are built 
out of symmetric and anti-symmetric combinations of identical Bogolubov factors. Thus 
the two kernels are intimately linked. For the case when the bath is a standard harmonic 
oscillator this interelationship can be written as a generalised fiuctuation-dissipation relation 



3 Evolution Operator and Master Equation 

In this section our goal is to calculate the evolution operator for the reduced density matrix 
and the master equation. The master equation is the evolution equation for the reduced 
density matrix. It provides a transparent means for sifting out the different physical processes 
caused by the bath on the system. First we must calculate the evolution operator in (2.7), 
which contains all the dynamical information about the open system. From this point on we 
shall put F{x) = X. 

The influence functional (2.10) and the corresponding influence action (2.8) can be written 
in a compact way 

v4[x, x'] = S[x] — S[x'] + 6A{x,x'), 

t t (3 1) 

SA[x,x'] = - 2 f ds f ds' A(s)/x(s,s')S(s') +i f ds ds'A{s)u{s,s')A{s') 

S[x]-S[x'] = ['ds{M{s)t{s)A{s) + -M{s)B{s)[^{s)A{s) + A{s)t{s)] 

2 (3.2) 

- M{s)Q\s)J:{s)A{s)} 

with the use of the 'center of mass' and 'relative' coordinates deflned earlier in (2.13). 

As pointed out by many authors 0, ^, ^, and in Sec. 2, the real and imaginary parts 
of A[x, x'] can be interpreted |^ as being responsible for dissipation and noise respectively. 
The imaginary part of (3.1) is determined by z^(s), the noise (or fluctuation) kernel. The 
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name becomes apparent when we realize that this term can be interpreted as coming from 
the interaction between the system and a stochastic force ^ that is hnearly coupled to the 
system and has a probability density given by Vl^] = exp{— ^(^z/)~^^}. On the other hand, 
the kernel n{s) in (3.1) is known as the dissipation kernel. The motivation for the name 
comes from the fact |p that it introduces a modification in the real saddle point trajectories 
of the path integral in (2.7). Strictly speaking only the non-symmetric part of the n kernel 
should be associated with dissipation. Thus, all the symmetric part can be absorbed in a 
non-local potential (that does not contribute to the mixing of the x and x' histories). There 
is no such symmetric part in the /i-kernel of our problem although it does appear in other 
cases 0. 



3.1 Evolution Operator 

The evolution operator given in equation (2.7) generates a non-Markovian dynamics since 
it fails in general to satisfy the relation 

for the reason that the operator Jr{t2iti) depends on the state of the system at time ti, 
unless that time is the one for which the system and the environment were decoupled. The 
non-Markovian behavior is, in fact, a direct consequence of the non-locality of the influence 
functional. 

Our task is to compute the evolution operator 







1 y A- f- 

1 '—^11 '^t 


)= DY. DA exp 





^^[S(.),A(.)] 



(3.3) 



Let us schematically describe how to compute the path integral in (3.3). We start by 
reparametrizing the paths, writing 

E(s) =x+(s) + S,/(s) 

(3.4) 

A{s) = x_{s) + A,i{s) 

where the "classical paths" f ^ ^ are solutions to the equations of motion derived from 



the real part of A[T,, A], and x± are the deviations from the classical paths. The equations 
governing these functions are 



B{s) M{s)B{sy 
~Y~ ^ 2M{s) 



M{s) 



ds'fi{s, s')S,z(s') = 



and 



Se/(t) = Sj 



(3.5) 
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and 



M{s) 



2 ■ 2M(s) J ' M{s) 
Aa{ti) = Ai, and Aci{t) = A/. 

After the path-reparametrization, (3.3) can be rewritten as 



where 



Jr{T^f,Af,t I Ei,Ai,ti) = Z(t,tj)exp 
Z(t,U) 



h 



A[E^{s),A^{s)] 



t■,x-^.=0 pt;x—=0 

Dx+ / Dx- exp 

ti;x-{-=0 Jti;x-=0 



-A[x+{s),x^{s)] 



1 /■* /■* 

We can write the classical solutions Tld and in terms of the elementary functions 



f'ds f ds'[x4s)Aa{s'Hs,s')] 



(3.6) 
(3.7) 

(3.8) 



Aa{s) = AiVi{s) + ^Ms) 
which satisfy the boundary conditions 

ui{s = ti) = 1 = U2{s = t) 
ui{s^t)^ = U2{s^ti) 

Vlis = ti) = 1 = V2(s = t) 

vi{s = t) = = V2{S = ti). 

Now setting 

bi(t, ti) = M(t)u2(t) + , 62(t, ti) = M(ti)u2(ti) 



b3{t,ti) = M{t)ui{t), hi{t,ti) = M{ti)ui{ti) + 
where the dot denotes the derivative with respect to s and 



7V/(t,)i?(t^) 



1 * * 

aij(t,ti) = / (is / (is't'i(s)i/(s,s')^j(s') 

i + Ojj J J 



ti ti 



we get 

Jr{xf,x'j:,t I a;j,a;^,tj) = Z{t,ti)exp 



h 



{h^fAf - b2^fAi + hJ^iAf - 64SiAi} 



X exp 



--{aiiA^ + auAiAf + assAj} 



(3.9.a) 
(3.9.6) 

(3.10.a) 
(3.10.6) 



(3.11) 



(3.12) 



(3.13) 
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The evolution operator (3.13) must preserve the normahsation of the density matrix. By 
requiring that Tr{p) = 1, (2.5) imphes 



We therefore find that 



(3.14) 



3.2 Master Equation 



We now proceed with the derivation of the master equation from the evolution operator 
(3.13) using the simplified method of Paz [Q. We first take the time derivative of both sides 
of (3.13), multiply both sides by p,.(Sj, Aj, tj) and integrate over Sj, Aj to obtain 



Pr(S/, A/,t) 



Pr(S/, Aj,t) 



+ -A/63 / dAidEiEiJrPr{Tji,Ai,ti) 



(ziaSj + diaA^) J dAd^,AiJrPri^^,\,U) (3-15) 



h 



- -hi I dAidT^iT^iAiJ-rPriT^i, Ai.ti 



On 

h 



J dAidEiA^iJrPr{^i,Ai,ti). 



Here the dot denotes derivative with respect to t. We can perform the integrals in (3.15) by 
using 

ihdJr ^ biAf 



62 <9Sj 



02 



^i3r — , 



dj 

h-^^+{Aiau + 2Afa22)Jr, 
Uh d 61 A A 

\h2Wf^~b^) 



X 



b3dAf bs' &3 / 



(3.16a) 
(3.166) 

(3.16c) 



The derivation of the master equation simplifies greatly with the use of the following relations 

(3.17) 



Ui{s) = Wi{s) - W2{S) 77^, U2{s) 



W2{t) 



W2{t) 
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In order to satisfy the boundary conditions, (3.10a), we require 

Wi{ti) = W2iti) = 1, W2i0) = M^) = ^■ 

In this representation we can show that 
b4 



'2^2, 



With these relations the master equation reduces to 



2M{t) ^dx^ dx'^) + y^^H^^ + ^ a^j 



+ 



M{t) 



nl^{t, ti) {X^ - X'^) + ih^] Prix, X', t) 



d d 

- ihT{t,ti){x - - ^) pr{x,x\t) 
+ iDpp{t,ti){x — x'Y Pr{x,x',t) 

- fi{D^p{t,ti) + Dp^{t,ti)Yx - x')(^-^ + ^) pr{x,x',t) 



92 



where 



bih h , B\t) b2B{t) 



r{t,ti 



f 0-12 _ On 



b2 \M{t) 



M{t)h M{t) 4 

2[b3 b2) 

ail \ 261 
b2 



262 



• ,0^3 , bib^ 

— — — a22 — Oj22 -r ^7—022 + '^127— [- 

M[t) 63 62O3 



. bi 

ai2T- 
02 



ai2 _ Q22 _ ^3^12 _'^( Q12 _ On 
62 M{t) bsb2 62 \M{t) 62 , 



' 62 \M{t) 62 , 



(3.18) 



(3.19) 



(3.20) 



(3.21) 
(3.22) 
(3.23) 
(3.24) 
(3.25) 



The dot in these equations denotes the derivative with respect to t. We can rewrite the 
master equation in the operator form which may be easier for physical interpretation. We 
find that it becomes 



d 

ih—pr{t) ^[Hren, p] + iDpp[x, [x, p\] + iD^xlP, [P, p\] 

+ iDxp[x, [p, p]] + iDpxlp, [x, p]] + T[x, {p, p}] 



(3.26) 
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where 

Hren = " ^{PS: + Xp) + ^fi,en(t)x'. (3.27) 

From the master equation we know that D^x and Dpp generate decoherence in p and x 
respectively and F gives dissipation. The master equation differs from Paper I by more 
than changing the kernels. The factor ai2/M{t) — aii/62 vanishes only when the dissipation 
kernel is stationary (i.e a function of s — s') and also when the system is a time independent 
harmonic oscillator. When this happens fi(s) = U2{t — s) and we have vi{t) = —b2/M{t). 
We see from (3.19) that the factor ai2/M{t) — aii/62 is zero in this case. All the diffusion 
coefficients contain this factor and D^x depends solely on it. Thus Dxx arises purely from 
non-stationarity in the dissipation kernel and system. 

The coefficients Dxx, Dpp, Dxp and Dpx promote diffusion in the variables p^, and xp+px 
respectively. This can be seen by going from the master equation to the Fokker-Planck 
equation for the Wigner function |l|, ^ . The Wigner function is defined by 



1 r°° , A A 

t) = ^ e^^^/^S - y + -)ciA. (3.28) 



where S, A are defined in (2.13). We can show that the Wigner distribution function from the 
master equation (3.26-7) (with B{t) = 0) obeys the following Fokker-Planck type equation 



d 

—Fw{^,p,t) 



M{t)dE ' 2 dp ' ^'^'dp'' ~ ' dp^ 

- hDxxit)—^ + 2(pxp{t) + Dpxit))^^^ Fw{T.,p, t). 



(3.29) 



4 Simple Examples 

4.1 Squeezed Thermal Bath of Static Harmonic Oscillators 

This is the simplest case treated before in Paper I and II. In this case the bath modes have 
time independent coupling constants with the Lagrangian 

m = \[e-uw (4.1) 

From (2.1) m„ = 1, 6„ = and tu^ = u"^. Substituting these into (2.21) and solving (2.20) 
(with n = uj) one obtains 

a = 6-'"^', /3 = (4.2) 
where a = 1 at the initial time t = 0. Substituting these into (2.18-19) one obtains 

POO 

n(s,s') = — / duj I (lu) sin uj(s — s') (4.3) 
Jo 
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and 



s') = / duj coth — — — I{uj) cosh2r(co') cos[co'(s — s')] 
Jo yikBi J L 

- smh.2r{uj)cos[2(j){uj) - uj{s + s) . (4.4) 

This is a simple generalisation of previous studies in that we have a squeezed thermal initial 
state rather than a thermal state. There are two distinct contributions to the noise 
kernel for an initially squeezed bath. The first term represents a change in the spectrum of 
the stationary vacuum noise. The second term has a new feature which is a non-stationary 
contribution to the noise kernel. This is expected since the fluctuations of a squeezed vacuum 
mode oscillate between conjugate observables. 

As (s + s') — >■ oo the second term in (4.4) tends to zero. Thus the nonstationarity in the 
noise kernel is a transient effect for the initial squeezed bath. For an initial squeezed bath 
with thermal spectrum the late time noise kernel would tend to that of the usual thermal 
state. This is because at late times , the noise kernel v loses track of the initial phase 
distribution 9{uj). This is, however, not true for the master equation diffusion coefficients. 
Equations (3.23-25) show that the diffusion coefficients depend on the noise kernel in a non- 
local way in time. It may be interesting to compare the timescales in which the semi-classical 
system and the full quantum system forget the 0(0;) initial condition in the bath. 

Although we have considered only single mode squeezed initial states our results can be 
easily extended to two- mode squeezed initial states [^]. This will change the noise kernel 
(4.4) but not the dissipation kernel (4.3) which remains independent of the initial state. 
Since the influence functional (2.10) is unchanged the exact forms for the evolution operator 
and master equations in Sec. 3 will stay. These results could then be used for an accurate 
description of systems coupled to an initially squeezed bath[^, |5^ . 

If we set the initial squeezing to zero we regain the results of Paper I. For completeness 
we will summarise the simple analytical results obtained previously. In this case the noise 
and dissipation kernels are functions only of s — s' . They can always be related by some 
integral equation known as the fluctuation-dissipation relation (FDR) fl]]: 



u{s) = J ds' K{s - s') 7(s') (4.5) 
—00 

where the kernel K{s) is 

K{s) = / — uj coth. -jShuj cos Lus (4.6) 
J n 2 



and fi{s) = -^lis). In the classical or high temperature limit, the kernel K is proportional to 
the delta function K{s) = 2kBT6{s) and the FDR is equivalent to the well known Einstein 
formula. 
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An interesting case is an environment which generates an ohmic spectral density 



= -IqMu. (4.7) 

TT 



With a discrete high frequency cutoff A, 

, , 2 , ^ (i /sin As 
/i s = -7oM— 

TT as \ s 

27oM^5(s), as A ^ oo. (4.8) 

as 

In this case for a harmonic oscillator system (2.17) becomes 

X{s) + 27oX + nlX = -e(t) (4.9) 

where Qr = ^ - -^o-^- We see that the ohmic environment is special in that it gives local 
dissipation in the infinite cutoff limit. 

Theoretically, the meaning of renormalization can be understood as follows []1|: We can 
rewrite the action as 

S= fds{lM{x'-Q'x')+Y.[lm^ql-lm.^ul{q + ^xr + l^x%} (4.10) 
Jo ^2 V 2 2 mnUJ^ 2m„u;^ J 

The last term can be viewed as a frequency counter term Q'^ arising from the interaction of 
the Brownian particle with the bath oscillators 

^^^ = -— EA = -/^-^- (4-11) 
The bare frequency is thus modified into a renormalized frequency fi^ given by 

^il = n^ + nl (4.12) 

Another interesting case is the high temperature limit. If we consider the temperature 
to be such that << A~^ and then let A — > cx) (the order in which we take the limits is 
important), the noise kernel (4.4) is simplified to 

z/(s) = 6{s). (4.13) 

In this case we see that the noise is white with an amplitude A'^oMhsT, and (4.9) reduces 
to the Nyquist formula. In the ohmic, high temperature and infinite cutoff limit the master 
equation coefficients can be calculated. Using (3.5) we find that, for a time independent 
harmonic oscillator system, ui and U2 must satisfy 

u{s) + 27om(s) + n^uis) = -47o5(s)m(0). (4.14) 
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The solutions satisfying the appropriate boundary conditions (with tj = 0) are 



ui{s) 



sin[fi(s - t)]e-^o" 
sin Clt 



sm[Qs]e' 



-7o(s-*) 



sin ilt 



where fl"^ = fl^ — 7q. Applying these to (3.11) we find 



hit) 
hit) 



sin Vtt sin Vtt 

-hi (t) = M(7o - cot nt) . 



(4.15) 

(4.16) 
(4.17) 



Since 64 is discontinuous before and after t = (due to the kick) we have taken the average. 

The results (4.16-17) are exact in the infinite cutoff limit of an ohmic environment. 
This is a local approximation which has been shown to be good for timescales greater than 
the inverse cutoff [^. Equations (4.16-17) depend only on the dissipation kernel which is 
unchanged by initial squeezing in the bath. Thus these equations can also be applied to 
more general situations. 

Using the noise kernel (4.13) and the fact that fi(s) = U2{t — s), ^2(5) = Ui{t — s) we 
can calculate aij and find that the master equation coefficients to be 

nrenit)=nr, r(t) = 7o, D,p{t) = D,,{t) = 0, Dpp{t) = . (4.18) 



For decoherence studies under these and other environmental conditions see 20 



4.2 Bath of Upside Down Oscillators 

This is the next simplest case. In this case the bath modes have the Lagrangian. 



(4.19) 



From (2.1) m„ = 1, 6^ = and 
(with K = uj) we obtain 



^. Substituting these into (2.21) and solving (2.20) 



Ciuj(t) = coshcijt, Puj{t) = —isinhut (4.20) 
where a = at t = which is our initial time. Substituting these into (2.18-19) we obtain 

H{s, s') = - dujl{uj) sinhcj(s - s') (4.21) 
Jo 

and 



z/ s, s 



^°°d..coth(^j 



cosh 2r(uj) coshct;(s -|- s') 



sinh2r(c(j) cos2(f){uj) coshct;(s — s' 
sinh 2r{uj) sin20(c(j) sinhco'(s -|- s') 



(4.22) 



This case can be used as an amplifier model in quantum optics and electronics pB . 



20 



5 Particle Detector in a Scalar Field Bath 



The formalism developed here can be used to study quantum statistical processes in cos- 
mological and black hole spacetimes. The model (2.1) can be used to depict a particle 
detector in motion, or an observer near a black hole. It can also be used to describe the 
non-equilibrium dynamics of homogeneous and inhomogeneous modes (density fluctuations) 
of the infiaton field or gravity wave perturbations (which in the linear approximation obey 
the wave equation of a massless, minimally coupled scalar field) in the early universe. 

In this section we will show how a general real scalar field in an expanding universe is 
reduced to a sum over quadratic time dependent Hamiltonians. The action for a free massive 
(m) scalar field in a curved spacetime with metric g^i, and scalar curvature R is given by 



S 



C{x)d^x 



-9 



(5.1) 



where \/y denotes covariant derivative, and is the field coupling constant (^^ = 0,1/6 
respectively for minimal and conformal coupling). In the spatially- fiat Robertson- Walker 
(RW) metric 



ds^ = a^{rii)[drf' — dxf] 



(5.2) 



we can write 



Cix) 



(5.3) 



where a dot denotes derivative taken with respect to conformal time r] = J dt/a. If we rescale 
the field variable x = (^^y this becomes 



C(x) 



- T.ix.r - (m'a' + {QU - 1)- - (1 - Qi 



d f a 



di] \a 



-X 



(5.4) 



where the last term is a surface term. [| 

If we confine the scalar field in a box of co-moving volume (fixed coordinate volume), 
we can expand it in normal modes 



X{x) 



cos k ■ X + Oj- sink ■ x\ 

k k ' 



(5.5) 



which leads to the Lagrangian 



.a\2 



2(1 - 66)7^1^1 - I ^ + « + (6^'^ - I 



(5.6) 



^The part of the surface term proportional to has been added in by hand. The surface term ensures 
that the second derivative of the scale factor doesn't appear in the Lagrangian |Q. This is necessary to 
have a consistent variational theory when the scale factor is treated dynamically rather than kinematically 
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where k — \k\ and Llrj) — J jC{x)(Px. Canonical momenta are 
Defining the canonical Hamiltonian the usual way we find 



PS = ^=«-(l-6e.)-^E. (5.7) 



2' 



Pf + (1 - 6e.)^(p|?| + <lM) + {k' + m-a- + 6e.(6e. - 1)-, qf 



(5.8) 

where the sum is over positive k only since we have an expansion over standing rather than 
travelling waves. 

The system is quantized by promoting {pZ, g2), (p'^^, qZ) to operators obeying the usual 
harmonic oscillator commutation relation. Thus the amplitude functions of the normal modes 
behave like time-dependent harmonic oscillators. (The Hamiltonian is not unique but is a 
result of our time coordinate and choice of canonical variables.) 

The above shows that a scalar field can be represented as a bath of parametric oscillators. 
In order to study the noise properties of the quantum field, we now introduce an interaction 
between the system, which can be a particle detector or a field mode, and the bath, the 
scalar field. 



5.1 Spectral Density of a Scalar Field 

Consider the general form of interaction between the system harmonic oscillator r, and a 
scalar field x oi the form 

^int{x) = -erxix)S{xo). (5.9) 

They are coupled at the spatial point Xq with coupling strength e. We want to derive the 
spectral density function for this field I{cu) — J2 ^{'^ —'^n)c\l2Kn. Integrating out the spatial 
variables we find that 

Lintir]) = j Cint{x)(fx = -erx(fo, n) (5-10) 

where 

x{xo,ri) = J ^Y.[qi cosk ■ xo + sink ■ xo]. (5.11) 

' k 

Comparing this with (2.1) we see that each set of modes has the effective coupling constants 

ci ^ '^l^e cos k ■ Xq, = ^y/^esinA; • fo- (5-12) 

In the continuous limit the oscillator label n is replaced by k. Adding the spectral densities 
from both the ± sets of modes we obtain 

k 
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where oj is replaced by k. In the continuous hmit: ~^ (^) / Writing d?k = 
k"^ sin 6dkd6d(f) and integrating between the hmits 0[27r, 0] and 6[k/2,0] (remembering we 
only include half of the modes) j^yi Io° k'^dk, we get 

For a two-dimensional scalar field we get 

I{k) 



e2 



2'KKk 

5.2 Accelerated Observer 

We consider a two dimensional massless scalar field $ in flat space with mode decomposition 

[J 

= Y ^ X![^a!' cos kx + sin kx]. (5.15) 

' k 

The Lagrangian for the field can be expressed as a sum of coupled oscillators with amplitudes 
for each mode 



o- k 

Now consider an observer undergoing constant acceleration a in this field with trajectory 

x(r) = - cosh ar, s(r) = - sinh ar. (5-17) 

We want to show via the influence functional method that the observer detects a thermal 
radiation. This effect was first proposed by Unruh |33], as inspired by the Hawking effect 



| 31| for black holes. Let us see what the spectral density is. The particle- field interaction is 

Cintix) = -er<l>{x)6{x{T)) (5.18) 

where they are coupled at the spatial point x(r) with coupling strength e and r is the 
detector's internal coordinate. Integrating out the spatial variables we find that 

L.ntir) = J Cint{x)dx = -er$(x(r)). (5.19) 

Comparing (5.19) with (2.1) we see that the accelerated observer is coupled to the field with 
effective coupling constants 

^'^^^^ ^ ey^cosA:x(r), c-(s) = e^sinA;x(r). (5.20) 
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From (2.26) the spectral density in the discrete case is given by 

6{k - KKirKir') 



+- 



(5.21) 



where we have to include the sum over both sets of modes and we have put Kn = = \kn\- 
This ensures that fn{si) = in (2.21). Making use of (5.20) and Y.n ^ ^ I dk we find that 
(5.21) becomes 

I{k, r, r) = I{k) cosk[x{T) - x{r')] (5.22) 

where I{k) = is the spectral density of the (2-dim) scalar field seen by an inertial 
detector. From (4.3) and (4.4) we can write, using an initial vacuum state, 



C(s(r), s(r')) = z/(s, s') + ifx{s, s') 
We can rewrite this as 

C(r,r') 



dk I{k, r, r')e 



/\ ~iuj[s{T)-s{T') 



(5.23) 



1 

2 Jo 
1 r'^ 



+ 9/ 

2 Jo 



dk /(A;)e~*'^[^*^'')~^*^'''^+''(''^~'^('^')1 
dk /(A;)e~*'''^'*^^'^"'''*^^^+''^^''~*^^'^1 



(5.24) 



dk' I{k') exp (-2ifc'e''^sinh[aA]/a 



which upon using (5.17) can be written as 

1 

^(^■^ ' - -2 1 

+ exp (^-2ik'e~''^ sinh[aA]/a 

where 2S = r + r' and A = r — r'. Making use of P| 

g-iQsmh(x72) ^ _ / c^;/f^2i,y(a;)[cosh(7rz/) cos(ra) — z sinh(7rz/) sin(z/a;)] 

TC Jo 

where Kn is a Bessel function of order n, we find that (5.25) becomes 

'irk 



C(r,r')= / dkGik) 



coth ( — 1 cos k{T — t') — i sin k{T — r') 



(5.25) 
(5.26) 

(5.27) 



where 



2 r / 

G{k) = —smh{7rk/a) dk' I{k') K2ik/a (^k'e"^ /c 
na Jo L ' V 

+ if2^fc/a (2A;'e-'^^/a)] = I{k). 



(5.2^ 



In deriving this we have used the integral identity 

dx x^'KJax) = 2^-^a-'^-ir 



(5.29) 

and the properties of gamma functions. Comparing (5.27) with (4.3-4) we see that a thermal 
spectrum is detected at temperature 

ksT = ^. (5.30) 
Zn 

This was first found by Unruh |^3[ and stated in this form recently by Anglin [pl|]. 
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5.3 Hawking Radiation in Black Holes 

Consider the metric of a two-dimensional uncharged black hole with mass m 



2m\ 
r 



r J 



In the Regge- Wheeler coordinates 

r* = r + 2m In |r/ (2m) — 1| 

this can be written as 



ds' 



2m 



{dt'-dr*'). 



The Kruskal coordinates are defined by 

-r* -t 



t — r* — —4m exp 
With this the metric becomes 

ds^^ 



4m 



, t + r* = 4m exp 



r* + t' 
4m . 



2m 



(5.31) 

(5.32) 
(5.33) 

(5.34) 
(5.35) 



Since the metric (5.35) is conformal to flat space, the fleld theory is equivalent to that of 
fiat space. Thus a detector at constant Kruskal position f* will have an influence functional 
identical in form to that of an inertial detector in flat two-dimensional spacetime in Kruskal 
coordinates. However we are interested in a detector at constant r*. In this case wc see from 
(5.34) that constant r* is effectively an accelerating detector in Kruskal coordinates since 

f*{t) = 4me''*/(^"^) cosh[V(4m)]. (5.36) 

We also want to express the influence functional in cosmic time t which from (5.34) is 

t{t) = 4me'^*/(^"*) sinh[V(4m)] (5.37) 

for the detector at constant r*. This case is now similar to the accelerating observer and as 
in Sec. 5.2 the spectral density is 



I{k,t,t') ^ I{k) cos k[f*{t) -f*{t')] 



(5.38) 



where I{k) — ^ and u — \k\. With this spectral density we can write for a massless scalar 
fleld in a two-dimensional black hole spacetime 



C{t,t')=iy{t,t')+iii{t,t') 



+ 



2 Jo 



2 Jo 
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dk /(A;)e"'''[''*W"*'*(*')+*(*)"*(*')l 
dk 7(A;)e~^'=[^*(*')~*'*(*)+*(*^~*(*')l. 



(5.39) 



Comparing (5.39) and (5.24) we see that this case is identical to the accelerated observer if 
we identify a = l/(4m). The factor involving r* can be absorbed into the definition of k. 
Hence we can rewrite (5.39) as 

({t,t')= dk I (k) [coth{A7rmk) COS k(t-t') -i sin k{t-t')]. (5.40) 
Jo 



Comparing (5.40) with (4.3-4) we see that a thermal spectrum is detected by an observer at 
constant r* at temperature 

kBT=^. (5.41) 



In the two dimensional case the detector response is independent of its position r*. This 
will not be the case in four dimensions. 



5.4 Hawking Radiation in de Sitter Space 



We now illustrate how the Gibbons-Hawking result can be obtained from the influ- 
ence functional method. These examples are also of practical use for describing the non- 
equilibrium dynamics of the homogenous and inhomogenous (density fluctuations) modes of 
the inflaton field in the early universe gT|, |2|, ||, ||, ||, 



5.4.1 Massless conformally coupled field 

Consider now a four-dimensional spatially-flat Robertson- Walker (RW) spacetime with met- 



ric 



(5.42) 



For this metric the Lagrangian density of a massless conformally coupled scalar field, defined 
by (5.1), is 



1 













a J 



(5.43) 



where a dot denotes a derivative with respect to t. Decomposing $ in standing wave normal 
modes we find (after adding a surface term) 



a? 



a 



(5.44) 



where k=\k\. If we wrote the Lagrangian in terms of conformal rather than cosmic time we 
see from (5.6) that we would have obtained a bath of stationary oscillators. Our kernels 
would then be (4.3) and (4.4) but written in conformal time. If we were to rewrite these 
kernels in cosmic time we would get the same kernels as those by starting with a Lagrangian 
written in cosmic time as we are doing here. 
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The detector-field interaction is of the same form as (5.9) (with x replacing $) and we 
find that with = k (5.14) gives 



I{k) 



e 

2^ 



k. 



Using the Lagrangian (5.44) we find from (2.20) that the Bogolubov coefficients are 



a 



2a 



-ikr) 



/3 



'1-a' 
2a 



-ikrj 



(5.45) 



(5.46) 



where r] = J^^dt/a(t) with a(tj) = 1. Using these we find that the noise and dissipation 
kernels (2.18-19) are, for an initial vacuum state 



a{t)a{t') 



dk I{k)e 



-ik{ri~ri') 



(5.47) 



We will now specialise to the de Sitter dynamics where, in the spatially-fiat RW coordinati- 
zation 0], the scale factor has the time- dependence 



a(t) = e 



Ht 



(5.4J 



In this case -q = -^e'^^ with U = 0. If we define A = t-t', 2S = t + 1' we find that (5.47) 
becomes 

dk I{k) exp 



C(t,t') = e-^^^ 
Using (5.26) we find that 

CM- 

where 

G{k) 



e-"^ sinh(i7A/2) 



H 



(5.49) 



dk G{k) 



coth i^j "^os k{t — t') — i sin k{t — t') 



(5.50) 



4sinh(7rfc/if) 



2 



dk' I{k')K2^k/H{2k'e-''''/H) 



(5.51) 



We have again used the integral identity (5.29) and the properties of gamma functions. 
Comparing (5.50) with (4.3-4) we see that a thermal spectrum is detected by an inertial 
observer in de Sitter space at temperature 



ksT = — . 
2n 



(5.52) 



Cornwall and Bruinsma |4g] who considered the evolution of low momentum modes of an 
infiaton field coupled to a thermal bath in a de Sitter background also derived the infiuence 
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functional for a conformally coupled scalar field in de Sitter space. The noise and dissipation 
kernels they found in their Eq. (3.31) differs from ours since they did not add a surface term 
proportional to ^d- As a result they got nonstationary kernels when written in conformal time. 
As we pointed out previously a surface term is needed to give a consistent variational 
theory when the scale factor is treated as a dynamical variable. In this case we see from 
(5.6) that in conformal time conformal coupling with a bath of ordinary stationary oscillators 
gives the usual stationary kernels. In cosmic time these kernels lead to (5.50) which is still 
stationary, but shows the expected Gibbons-Hawking temperature. 



5.4.2 Massless minimally coupled field 

From (5.6) the Lagrangian for a minimally coupled massless field in de Sitter space is 



1 



(5.53) 



Solving (2.20) (with k„ = k) we find that the Bogolubov coefficients are 



1 - 



-ikri 



2kr] 



2kr] 



-ikrj 



Substituting these into (2.18-19) we find that 



((??, v') = J^ir], v) + if^iv, v') = dk I{k)e 

Jo 



-ik{v~v') 1 1 + k'^m' + - i) 

\ k'^rfr]' 



(5.54) 



(5.55) 



where I{k) is given by (5.45). We want to write this in terms of cosmic time given by 



r] = —jfC Following a similar procedure as before, we find 



C{t,t') = / dk G{k) 
Jo 



coth f j k{t — t') — i sin k{t — t') 



(5.56) 



where 



G{k) = I{k) 



' ^.H . , (Hit-t')\ , (ixk' 



(5.57) 



and we have ignored a factor e^(*+*') which gets cancelled by changing the integration measure 
from ?7 to if: in the infiuence functional. 

The imaginary part of (5.57) generates a contribution to the dissipation kernel of the 
form 



(5.5^ 



Inserting this into the infiuence functional (2.10) we see that it leads to a vanishing con- 
tribution to the infiuence functional. Similarly the imaginary part of (5.57) generates a 
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contribution to the noise kernel of the form 



2He^ 
2He^ 



sinh 



'H{t-t') 



J dk tanh (^'Jj^ smk(t — t') 



— sinh 



'H{t-t')\ cosA(t-t') 



t-t' 



H 



(5.59) 



where we have first integrated by parts and then used a standard integral. The first term 
in (5.59) will generate a vanishing contribution to the influence functional (2.10) since it 
involves an integral over a term oscillating infinitly fast. The second term in (5.59) can also 
be ignored since it generates only a zero frequency contribution to the noise spectrum. Thus 
the imaginary part of (5.57) can be ignored leaving a thermal influence functional at the 
Gibbons-Hawking temperature but with an effective spectral density of the form 



G{k) = I{k) 



1 + 



(5.60) 



We see in this spectral density the well known infrared divergence associated with massless, 
minimally coupled fields in de Sitter spacetime. 

Habib and Kandrup claimed |5^ , from a classical analysis, that a fluctuation-dissipation 
relation (FDR) would increasingly fail to hold as the physical period of oscillation increased 
over the expansion timescale of the universe. We suspect that the definition of FDR and 
its applicability in their work is more restricted than ours. We see that in both of these 
examples here the FDR (4.5) is exact despite the fact that the physical period of oscillation 
can be arbitrarily greater than the expansion timescale. This is consistent with the view of 
||16| , |l], that the FDR is a categorical relation as it is a description of the full backreaction 
of the environment on the system. 



6 Discussion 

Many physical problems can be modeled by a quantum Brownian particle in a parametric 
oscillator bath. We mention quantum optics, quantum and semiclassical cosmology and 
gravity. This paper aims to accomplish two goals: 

I. To derive the influence functional of a parametric oscillator bath, the evolution operator 
and the master equation for the reduced density matrix for explicit use in these problems. 

II. To relate the quantum mechanics of oscillators to quantum fields, thus providing a bridge 
from quantum stochastical mechanics to quantum field theory. This connection can benefit 
the former with the well-developed techniques of field theory (e.g., use of Feynman diagrams 
1^, ^) and enrich the latter with imparting a statistical mechanics meanings to many 



quantum effects R 11, 16 . 
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Two issues are discussed in this paper: 

A. The nature and origin of noise and dissipation in quantum fields 

B. The statistical mechanical meaning of quantum processes in the early universe and black 
holes. 

On the first issue we have discussed these problems: 

1) How to extract the statistical information of a quantum field. We couple a particle detector 
to the oscillator bath and study the detector's response to the fiuctuations of the field. We 
found that the characteristics of quantum noise vary with the nature of the field, the type of 
coupling between the field and the background spacetime, and the time-dependence of the 
scale factor of the universe. 

2) How to relate noise to particle creation. Parametric amplification of vacuum fiuctuations 
and backscattering of waves in the second-quantized formulation give rise to particle creation. 
By writing the infiuence functional in terms of the Bogolubov coefficients which determine 
the amount of particles produced, one can identify the origin of noise in this system to 
particle creation 13^, |35|, |36[ . 

On the second issue, we have studied the problem of 

3) quantum noise and thermal radiance. We illustrate how a uniformly accelerating detector 
in Minkowski space, a static detector outside a black hole and a comoving observer in a de 
Sitter universe observes a thermal spectrum. The viewpoint of quantum open systems and 
the method of infiuence functionals can, in our opinion, lead to a deeper understanding of 
black hole thermodynamics and quantum processes in the early universe 0. 



As further studies, the results obtained here can be useful for the following problems: 
a) Decoherence. The transition of the system from quantum to classical requires the dimin- 
ishing of coherence in the wave function. The noise kernel is found to be primarily responsible 
for this decoherence process. Decoherence can be studied by analyzing the magnitude of the 
diffusion coefficients in the master equation. The new result obtained here is useful for the 
analysis of decoherence where parametric excitation is present in the environment. This is 
the case when considering the quantum to classical transition of the wavefunction of the 
universe TB|, homogenous and inhomogenous modes (density fluctuations) of an infla- 
ton field EB, BTl |29[ or the primordial gravitational radiation background. For the case of 



density fiuctuations we can expect decoherence, dissipation and diffusion to have important 
consequences for the amplitude and spectrum of density perturbations. The relation of par- 
ticle creation and decoherence was one of the original physical motivations for this work. 
Indeed one of us has speculated [11| that in the early universe, vacuum particle creation and 
decoherence can be important at the same scale near the Planck time. We will address these 
issues at a later time. 

b) Backreaction The backreaction of these quantum field processes manifests as dissipation 
effect, which is described by the dissipation kernel pO|. In p6|, p5| we outline a program for 



studying the backreaction of particle creation in semiclassical cosmology in the open system 
framework. We use a model where the quantum Brownian particle and the oscillator bath are 
coupled parametrically. The field parameters change in time through the time-dependence of 
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the scale factor of the universe, which is governed by the semiclassical Einstein equation. We 
can derive an expression for the influence functional in terms of the Bogolubov coefficients 
as a function of the scale factor. The equation of motion becomes an Einstien-Langevin 
equation, from which a new, extended theory of semiclassical gravity is obtained. This, in 
our opinion, is necessary for furthering the investigation of quantum and statistical processes 
in curved spacetimes. We are currently pursuing these investigations from this viewpoint, 
c) A fluctuation- dissipation relation for non-equilibrium quantum fields. Sciama |^ first 
suggested that the thermal radiance in a uniformly accelerated observer (Unruh effect) and 
in black holes (Hawking effect) can be understood in terms of a fluctuation-dissipation rela- 
tion. This relation was also later derived for de Sitter spacetime via linear response theory 
by Mottola . These familiar cases all deal with spacetimes with event horizons and ther- 
mal particle creation. From earlier particle creation- backreaction studies in semiclassical 
gravity a general FDR was conjectured by one of us [16] for quantum fields in curved 
spacetimes. It corresponds to a non-equilibrium generalization of Hawking-Unruh effect to 
general dynamical spacetimes without event horizons. Such a relation can in principle be 
identified from the results of this paper. The interpretation of backreaction processes in 
terms of fluctuation-dissipation relations will be explored further in [Bl, RH . 
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A Influence Functional 



Here we describe the calculation of the influence functional. From (2.9) the influence func- 
tional is 

T{x, x'\ = Tr{ij{t, U)p,{ti)U'\t, ti)] (AA) 
where U{t) and U'{t) are the quantum propagators for the actions 

SE[q\+Sint[x{s), q] and SE[q\+Sint[x'{s), q], and x{s) and x'{s) are treated as time dependent 
classical forcing terms. 

Our first task is to determine the propagator for the action 

/ r 1 

SE[q] + Sint[x{s), q]^ J ds ^{^mn{s)(ql + hn{s)qnqn - '^n(«)?n)} 

ti L n 

+ Y.{-^^n{s)F{x{s))qn - C2n{s)F{x{s))qn (^-2) 
n 

- C3n(s)F{x{s))qn{s) - C4„(s)F(i;(s))g„) 

This interaction is the most general interaction possible which is linear in the bath. Dropping 
the n subscript the Lagrangian for a mode takes the form 



L{t) = lm{t) {e + Kt)qq - u;'{t)q') - q[c,{t)F{x{t)) + c,{t)F{xm 

-q[cs{t)F{x{t)) + c^{t)F{x{t))]. 
Defining the canonical momenta the usual way we find that 

Pc = = m{t)q + m{t)b{t)^ - [cs{t)F{x{t)) + c^{t)F{x{t))]. 

The Hamiltonian, H{t) = Pcq — L{t) then takes the form 



(^.3) 



{AA) 



Hit) 



2m(t) 
+ 



b(t), , m(t) ( 2/ X b'^{t)\ 2 

-^(Pc? + qPc) + -Y + ^ j 5 

h{t) 



Ci{t)F{x(t)) + C2{t)F{x{t)) - -^{csit)F{x{t)) + C4(t)F(x(t))) 
[cs{t)F{x{t)) + c,{t)F{xm ^ , [cs{t)F{ x{t)) + c,{t)F{x{t))Y 

' / . \ Pc ~^ 



{A.5) 



m(t) 2m{t) 
The system is quantized by promoting q, pc to operators obeying [q, pc] — ih. Then writing 



{A.Q) 
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we find that (A. 5) becomes 

H{t) = f{t)A + f*{t)A^ + h{t)B + d{t)a + d*{t)a^ + g{t) 
where A and B are defined in (B.2) and 



h{t) 



h ( K m{t)uj^{t) m{t)b^{t) 



d{t) 



2 \m{t) 

Ci{t)F{x{t))+C2{t)F{x{t)) 



+ 



+ 



hit) K 
I- 



9{t) 



2 ~m{t) 

c,{t)F{x{t)) + c,{t)F{x{t))f 
2m{t) 



Ak 



c^{t)F{x{t))+c^{t)F{x{t))] 



(A.7) 

(A.8) 

(A.9) 

(A.IO) 
(All) 



In appendix B we have derived the evolutionary operator generated by the Hamiltonian 
of (A. 7). It has the form 



f/(t,t,) = S(r,0)i?(^)D(p)exp 



— J g{s)ds+ I ds I ds'p{s)p*{s 



ti J t 



(A12) 



From the first two equations of (B.13) we see that the squeeze and rotation operators do 
not depend on x. Thus S = S' , R = R'. Using this fact, the unitiary nature of the operators 
in the propagator, the cyclic trace rule and the identity 



D{p)D{p') = D{p + p') exp 
we find that (A.l) becomes 



-{pp'* -p*p') 



(A13) 



T[x, x'] =Tr[pb{ti)D{p - p')] exp 



{pp'* — p*p — pp* — p'p' 



X exp 



ds / ds' 



p{s)p*{s')+p'*is)p'{s') 



h 



ds[g{s)-g'{s)] 



(A.U) 



Making use of the integral identity 

g{t)dt f h{t)dt 



a rt 



b Jb 



[g{t)h{t') + g{t')h{t)]dt'dt 



(A15) 



its possible to write 
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J^[x,x'] =Tr[/36(t,)L'(p-p')]exp[-- / ds[g{s) - g' {s)] 



X exp 



lf^dsf^ds\\p{s)~p\s)W{s')+p*{s')] (A16) 



+ W)+p\s'W*{s)-p*{s)]) 



We will now evaluate the influence functional for a squeezed thermal initial state. Our 
first task is to compute the trace in (A. 16). Our initial state is of the form 



Pb(ti) = S{r,(t))pthS\r,(t)) 



(A17) 



where pth is a thermal density matrix of temperature T defined by the thermal density matrix 
takes the form 



Pth 



1 — exp 



II exp 



-nhuo 



\n) {n\ 



(A18) 



and S'(r, 0) is a squeeze operator defined in (B.12). 
The trace in (A. 16) becomes 



Tr[f)i,{ti)t){p - p')] = Tr[pthS\r,(f))D{p - p')S{r,> 



Making use of |55 



S\r, (j))D{p)S{r, (p) = Z)(j9coshr + p* sinh 



re 



2i<t>\ 



equation (B.16) and ||70| 

Tr[pth exp(ta^ + ud)] = exp 

we find that 



tu , 
— coth 
2 



Tr[pb{ti)D{p-p')] = exp 



1 , f hu ' 
— coth — — — ■ 

2 \2kBT^ 



I {p ~ p') cosh r + {p — p'Y sinh re 



2i<f>\2 



(A19) 
{A.20) 

{A.21) 
(A22) 



Making use of the integral identity (A. 15) we can write 



Tr[pb(ti)D{p - p')] = exp 



coth f—^ I [ ds [ ds' 

2 \2kBT Ju Ju 



X |[p(s) - p' {s)\[p{s') - ]3'(s')]* cosh2r 

+ [p{s) - p\s)]*[p{s') -]3'(s')] cosh2r 
+ sm\i2re-'^'^[p{s) - p\s)][p{s') - p\s')] 

+ smh2re^''^\p{s)-p\s)Y\p{s') ~p\s')X 



(v4.23) 
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Now put into (B.17) 
where from (A. 10) 



d = uF{x) + vF{x) 



u 



and further define 




,C3 . C3 
Ci — IK — 

2 m 



, C4 . C4 

C2 - 0— - IK — 

2 m 



{A.2A) 



(A.25) 



(A.26) 



Using (B.17,A.23) and (A. 16) we find that the influence functional for a mode n takes the 
form 



Fn[x,x'] = exp 



where 



and 



2i 



t s 

ds J ds'[A{s)l^,n{s,s')^{s')+A{s)l^2n{s,s')^{s') 



+ A{s)l^3n{s, S')t{s') + A(S)/X4„(S, s')t{s') 
+ A{s>3nis, S')A{S') + A(5)z/4n(5, s')A{s') 



ds[g(s) - g'(s)] 



A{s) = [F{xis)) - F{x'is))], 2S(s') = [F{x{s')) + F(x'(s'))] 
A{s) = [F{x{s)) - F{x'{s))], 2t{s') - [F{x{s')) + F{x'{s'))] 

^i^nis,s') = ^[Uis)U*{s')-U*{s)U{s')] 

f^2n{s,s')^^[V{s)U*{s')-V*{s)U{s')] 

Hsn{s,s') = ^[U{s)V*{s')-U*{s)V{s')] 
nUs,s') = ^[V{s)V*{s')-V*{s)V{s')] 



J^ln{s,s') =— COth 



cosh2r{U{s)U*{s') + U*{s)U{s')) 



sinh2re-2'^t/(s)t/(s') - sinh2re=^^^t/*(s)t/*(s') 



{A.21) 



(A.28) 



(A.29) 
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and 



cosh2r{U{s)V*{s') + U*{s)V{s')) 



- smh2re-^'^U{s)V{s') - smh2re^'^U*{s)V*{s') 



1 



-3n(.,.')=^COth ^^^^ 



cosh2r{V {s)U*{s') + V*{s)U{s')) 



- sinh2re-2^'^y(s)f/(s') - smh2re'^''''V*{s)U*{s') 



Z/4„( S, S 



1 , f huj " 

coth — — — 
2h \2kBT 



cosh2r(F(s)\/*(s') + V*{s)V{s')) 



sinh2re-2^'^1/(s)y(s') - sinh2re=^*'^\/*(s)\/*(s') 



c^{s)F{x{s)) + Ci{s)F{x{s)) 



l2 



2m(s) 



(A.30) 



(A.31^ 



(A.32) 



(A33) 



Note that the g term in the influence functional can be absorbed into the system Lagrangian. 
Of course the total influence functional is an infinite product of influence functionals over n. 
Therefore if we define the spectral densitry as 



2 /t^t- 



we obtain the results of (2.18-19) where we have put C2, C3, C4 = 0. 
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B Propagator 

Consider the Hamiltonian 

H{t) = f{t)A + f*{t)A^ + h{t)B + d{t)a + d*{t)a^ + g{t) (5.1) 

where 

i=y, A^ = ^, B = d^a + l/2. {B.2) 

and [a, d^] = 1. We want to find the propagator for this general time dependent system. We 
make the ansatz 

(Jit^ U) = e-W^e^W^e^W^^e«(*)«ePW'i^e'-W. (S.3) 



It has proved by Fernandez that this is global. It must satisfy the evolution equation 
for the propagator 

H{t)U{t,U)=th^U{t,U) (5.4) 
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subject to the initial condition (j {ti, U) — 1. We find that the operators A, A', B, a, a) satisfy 
the following commutation relations 

[i, A^]^B = B\ [i, B] = i, [A\ B] = -2A^ 

[a, A'^] = a\ [a\A\ = -a 

[a,S]=a, [a\B]^-a} {B.5) 
[a, A] = [a^i^] = 0. 
Making use of the commutation relations and the operator relation 

,2 



g«Opg-«o _ p ^ ^[o, P] + —[O, [O, P]] + ... {B.6) 



we find 



qa 



{B.7) 



ei^a^ = (a^ + q)e 
e'^'d = {d- dh)e'^' 

e^^a = e-^ae^-^ 

e^^i^ = (it + + |/2i)e2'^ 
Substituting (B.3) into (B.4) and using (B.7) we find 

h = ih{x + zy) 

d = ih{q{l — yz)e^'^ + pye~^) 
d* = ih{pe^ - qze^) 
g ^ih{pq + r). 

Since the first three equations of (B.8) are independent of d and g the first three terms in 
the propagator (B.3) are independent of the last three. As it stands (B.3) is not necessarily 
unitary. Thus x, y, z must satisfy some further restrictions. If we write 

X = Ina, y = —^a, z = 13* /a {B.9) 



B.8 
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where 

a = e~*^ cosh r, P = —e~'^^'^ sinh r 
then we can write (B.3) as (using relations in [p5|] ) 

U{t,ti) = S{r,(j))R{e)e''^eP^\' 

where 2(j) = 2ip — 6 and 

R{e) = e-^'^, S{r, (j)) = exp[r(ie-2'^ - A^e^''^')]. 



(B.IO) 
(5.11) 

(5.12) 

S and R are called squeeze and rotation operators respectively ||55[. They are are both 
unitary as is required. Substituting (B.9) into (B.8) we find 

ha = —if* [3 — iha 

hp = ihp + if a 

hp = -i{dp* + d*a*) (S.13) 

q = -p* 
hr = —ig — hpq = —ig + hpp*. 

The first 2 equations of (B.13) completely determine a and /3. The last three determine 
p,q,r. Making use of 



IF,G] 



. " = e-e 2 (5-14) 
where F and G are any operators that satisfy [F, G] =constant, we find that (B.ll) becomes 

(5.15) 



where 
and 



U{t,U) = S{r,^)R{9)D{p)e-PP''/^e' 
D{p) = exp[pd^ —p*o] 



p{t,U) = -- / dt[d{t)p*{t) + d*{t)a*it)] 
h Jti 

r{t,ti) = f g{t)dt+ f p{t)p*{t,U)dt. 
a Jti Jti 



If we define 



P{t) = p^{t) + ip2{t) 

and use the identity, (A. 15), we find that (B.15) becomes 



(5.16) 

(5.17) 
(5.18) 

(5.19) 



lJ{t,U) = S{r,(t>)R{e)D{p)e^V 



i I ds I ds'[p2{s)pi{s') - pi{s)p2{s')\ 



X exp 



h 



g{s)ds 



(5.20) 



This form shows explicitly that the propagator is unitary. 
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